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This two-part contribution to tlie Proceedings of tiie Eightli Canadian Conference on General 
Relativity and Relativistic Astrophysics is devoted to the evolution of a massless scalar field in two 
black-hole spacetimes which are not asymptotically flat. 

In Part I (authored by Eric Poisson) we consider the evolution of a scalar field propagating in 
Schwarzschild-de Sitter spacetime. The spacetime possesses a cosmological horizon in addition to 
the usual event horizon. The presence of this new horizon affects the late-time evolution of the 
scalar field. 

In part II (authored by William G. Laarakkers) we consider the evolution of a scalar field propagating 
in Schwarzschild-Einstein-de Sitter spacetime. The spacetime has two distinct regions: an inner 
black-hole region and an outer cosmological region. Early on in the evolution, the field behaves as if 
it were in pure Schwarzschild spacetime. Later, the field learns of the existence of the cosmological 
region and alters its behaviour. 



Part I 
by Eric Poisson 

Introduction 

A theorem that establishes the uniqueness of the 
Schwarzschild black hole as the endpoint of gravita- 
tional collapse without rotation was proved by Werner 
Israel more than 30 years ago [Q, and the mechanism 
by which the gravitational field eventually relaxes to the 
Schwarzschild form was elucidated by Richard Price more 
than 25 years ago Given the venerable age of this 
topic, it is surprising that more can be said about it to- 
day. Yet, many papers on radiative falloff have been writ- 
ten in the last few years [|-|15). Most of the new develop- 
ments are concerned with rotating collapse, and how the 
gravitational field eventually relaxes to the Kerr form. 
The question we pursue in this two-part contribution is 
different. Focusing our attention on nonrotating black 
holes, we ask: How do the conditions far away from the 
black hole affect the relaxation process? In Part I we con- 
sider a black hole immersed in an inflationary universe. 
(This was first done by Brady et al. ||l^, and additional 
details can be found in Ref. (l^.) In Part II, William 
G. Laarakkers will consider a black hole immersed in a 
spatially-flat, dust-filled universe. 



Radiative falloff in Schwarzschild spacetime 

Price's result [|| can be summarized as follows. As 
a nonspherical star undergoes gravitational collapse, the 
gravitational field becomes highly dynamical, and the es- 
caping radiation interacts with the spacetime curvature 
surrounding the star. At late times, well after the initial 
burst of radiation was emitted, the gravitational field re- 
laxes to a pure spherical state. If 5g schematically repre- 
sents the deviation of the metric from the Schwarzschild 



form, then 5g ^ ^-(2^+2)^ where £ is the multipole or- 
der of the perturbation; the dominant contribution to dg 
comes from the quadrupole {£ = 2) mode. 

The inverse power-law decay applies to many other 
situations involving radiation interacting with the cur- 
vature created by a massive object. The simplest model 
problem which exhibits this behaviour involves a massless 
scalar field in Schwarzschild spacetime. In this context, 
the background geometry is not affected by the field $, 
which satisfies the wave equation 



[9 



(1) 



where gap is the spacetime metric, R the Ricci scalar 
(which vanishes for Schwarzschild spacetime), and ^ a 
coupling constant. Because the spacetime is spherically 
symmetric, the field can be decomposed according to 



(2) 



This leads to a decoupled equation for each wave function 
^i, and we can focus on a single mode at a time. 

The problem is formulated as follows. A pulse of scalar 
radiation (described by V'f) impinges on the black hole 
and interacts with the spacetime curvature, which creates 
a potential barrier fairly well localized near r — 2>M . The 
wave pulse is partially reflected and transmitted, and at 
late times, a tail remains. At such times, the field falls 
off as ^ ^-(2f+3)_ fpj^jg ig Price's power-law decay, and 
this behaviour is displayed in Fig. 1. 

A number of analytical and numerical studies of radia- 
tive dynamics l^-Q] have revealed that the inverse power- 
law behaviour is not sensitive to the presence of an event 
horizon. In fact, power-law tails are a weak-curvature 
phenomenon, and it is the asymptotic structure of the 
spacetime at radii r ^ 2M which dictates how the field 
behaves at times t ^ 2M. It is this observation that 
motivated our work: How is the field's evolution affected 
if the conditions at infinity are altered? 
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FIG. 1. Absolute value of the wave function tpe{t,r) as a 
function of time t, evaluated at r = 10 in Schwarzschild space- 
time. We use units such that 2M = 1. The cases £ = 0,1,2 
are considered, and the wave functions are plotted on a log-log 
scale. In such a plot, a straight line indicates power-law be- 
haviour, and a change of sign in the wave function is repre- 
sented by a deep trough. We see that the field's early be- 
haviour is oscillatory, but that it eventually decays according 
to an inverse power law. 

Radiative falloff in Schwarzschild-de Sitter spacetime 

To provide an answer to this question, wc remove the 
black hole from its underlying flat spacetime and place it 
in de Sitter spacetime, which describes an exponentially 
expanding universe. The Schwarzschild-de Sitter (SdS) 
spacetime has a metric given by 

ds^ = -f dt^ + f-^ dr^ + dVl^ , 

(3) 

/ = 1 - 2M/r - r^ja^. 

Here, a} = 3/ A, where A is the cosmological constant. 
(The SdS metric is a solution to the modified vacuum 
field equations. Gap + ^gap — 0, which imply R — 4A ~ 
12/a^.) The spacetime possesses an event horizon at r = 
Te ~ 2M and a cosmological horizon aX r = rc — a. We 
assume that re ^ rc, so that the two length scales are 
cleanly separated. 

We examine the time evolution of a scalar field in SdS 
spacetime; the field is still governed by Eq. (|l|), and it 
still admits the decomposition of Eq. (j^) . Figure 2 pro- 
vides a comparison between the behaviour of ■(/'£ in the 
two spacetimes (Schwarzschild and SdS). We see that at 
early times, the wave functions behave identically; the 
field has not yet become aware of the different condi- 
tions at r ^ rg. At later times, however, deviations be- 
come apparent. For i — Q, the Schwarzschild behaviour 
"00 ^ is replaced by the wave function changing sign 
dX t ^ 260, and settling down to a constant value at late 
times. For f = 1, the Schwarzschild behaviour ijji ^ t^^ 
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FIG. 2. Absolute value of the wave function tpi{t,r) as a 
function of time t, evaluated at r = 10 in Schwarzschild space- 
time (re — 1) and SdS spacetime (r^ — 1 and Vc — 2000). The 
cases ^ = 0, 1 are considered, and the wave functions are plot- 
ted on a log-log scale. In both cases, ^ = 0. 

is replaced by a faster decay which eventually becomes 
exponential. 

The field's exponential decay is confirmed by monitor- 
ing its evolution up to times t > r^. If ^ = 0, we find 
that 11)1 ^ 5"^"'=* at late times where Kc — 1/rc is 
the surface gravity of the cosmological horizon. 

A rich spectrum of late-time behaviours is revealed 
when ^, the curvature-coupling constant, is allowed to be 
nonzero. Figure 3 shows the time-dependence of ipo for 
several values of For ^ smaller than a critical value ^c, 
the field decays monotonically with a decay constant that 
increases with increasing ^. When ^ > ^c, however, the 
wave function oscillates with a decaying amplitude. As ^ 
is increased away from the critical value ^c, the frequency 
of the oscillations increases, but the decay constant stays 
the same. 

This qualitative change of behaviour as ^ goes through 
is quite remarkable. It can be explained with a detailed 
analytical calculation that will not be presented here (see 
Ref. ||l^). This calculation reveals that at late times, the 
field behaves as ^pi ~ g-pi^d^ where 

P='+|4v^+0(li). ,4) 

This relation implies that p becomes complex, and ipi 
oscillatory, when ^>^c = 3/16. 

Part II 
by William G. Laarakkers 

The spacetime 

The background spacetime in which the scalar field's 
evolution is followed is the Schwarzschild- Einstein-de Sit- 
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FIG. 3. Absolute value of the wave function tpo{t,r) as a 
function of time t, evaluated at r = 10 in SdS spacetime 
(re = 1 and = 100). Several values of ^ are considered, in 
the interval between ^ = and ^ = | • The wave functions are 
plotted on a semi- log scale, in which a straight line indicates 
exponential behaviour. 



ter spacetime. Qualitatively, it can be described as fol- 
lows. The idea is to start out with a spatially- flat, ex- 
panding, dust-filled universe. Then a ball of dust is 
"scooped" out, which leaves behind a spherical vacuum 
region. The dust that was removed is replaced by a 
Schwarzschild black hole, which is placed in the mid- 
dle of the vacuous region. This produces a spacetime 
with two distinct regions. The inner (black hole) re- 
gion is described by the Schwarzschild metric, and the 
outer (cosmological) region is described by the Friedman- 
Robertson- Walker (FRW) metric (see Fig. ||) . 

There are two important things to note about this 
spacetime. First, if the mass of the black hole is the 
same as the mass of the dust that was scooped out, the 
metric will be smooth across the boundary separating 
the two regions of the spacetime. Also, since the dust is 
pressureless it will not flow across the boundary, and the 
boundary itself will be co-moving with the universe. 

Because the specific finite-difference equation used in 
the numerical work requires the use of null coordinates 
(see ||]), the metrics of the two regions must be put in 
double-null coordinate form. For the black hole region 
the metric is written as 
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2M 

r 



dudv + dfl^ . 



(5) 



Here, u and v are ingoing and outgoing null coordinates, 
and r is defined implicitly by r -I- 2Mln(r/2A/ — 1) = 
(v — u)/2. In the cosmological region the metric takes 
the form 



a^{u*,v*) {-du*dv* + dn^), 



(6) 



where u* and v* are ingoing and outgoing null coordi- 
nates of the FRW spacetime, different from u and v. The 
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FIG. 4. Schematic of the Schwarzschild-Einstein-de Sitter 
spacetime. Note that the boundary between the two regions 
is expanding outwards, co-moving with the universe. 

FRW radial coordinate is x — \ — u*) . The scale fac- 
tor a is given by a{u*,v*) = jqC{u* where C is a 
constant that depends on the mass M of the black hole 
and the density of the dust. 

The first task is to find one coordinate system that 
can describe both regions of the spacetime. This is re- 
quired so that a single wave equation valid over the entire 
spacetime can be constructed. Since it is known that the 
metric is continuous across the boundary we can evaluate 
the metric induced on both sides of the boundary hyper- 
surface, and set them equal. This construction allows us 
to find the ingoing Schwarzschild coordinate it as a func- 
tion of the ingoing cosmological coordinate u* , and the 
outgoing coordinate w as a function of v* . Thus we now 
have a single coordinate system covering both regions of 
the spacetime. 



The wave equation 

The wave equation that governs the evolution of the 
scalar field is one without curvature coupling (equivalent 
to setting ^ = (see Part I and [^). Thus the massless 
scalar field <i> obeys the equation 



0. 



(7) 



The spherical symmetry of the problem allows us to de- 
compose the field in terms of spherical harmonics, and 
then to evolve only the part of the field that depends on 
the null coordinates. Thus the field can be decomposed 
as 



im 



(8) 



where R — r, -tpe — 'ipi(u,v) in the black hole region, 
and R ~ ax, i^i — iptiu*, v*) in the cosmological region. 
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When all quantities are expressed in the starred coordi- 
nate system, each wavefunction ^pi satisfies the equation 



du*dv* 



(9) 



Comparison of \|/ at Horizon for 1=0 yiQ=^ 



where the potential V takes a different form depending 
on which region of the spacetime the field lies: 



du dv 

V^Schiid — -rT-j-TJ 
du* dv 



£{£ + !) 



2M 

r 



VpRW — 



4e{i + i) 



(10a) 
(106) 



Results 

The numerical code evaluates the field on the event 
horizon, on the boundary between the two regions of the 
spacetime, and at future null infinity. Our discussion 
here will be restricted to the value of the field on the 
event horizon for the £ — and £ = 1 modes. The evo- 
lution was started at a "late" time, meaning that the 
boundary has expanded far enough that it can be clearly 
seen that the field initially behaves as it would in pure 
Schwarzschild spacetime. For both modes considered we 
see in figures || and || that the field first exhibits quasi- 
normal ringing followed by the well known power law de- 
cay (see, among others, [^). However, at a certain time 
in the evolution, the field's behaviour deviates from the 
behaviour exhibited in the pure Schwarzschild case. The 
point at which the field changes behaviour corresponds 
to the time at which information about the existence of 
the cosmological region reaches the event horizon. 

As the wave packet falls towards the event horizon (ap- 
proximated by u* = Mmax, whcrc Mmax IS the largest value 
of u* in the numerical grid) it encounters the localized 
potential (dashed line — see Fig. 0). Part of the wave 
is transmitted through the barrier and reaches the event 
horizon, and part of the wave is back-scattered by the 
potential. The reflected wave heads out towards the cos- 
mological region (to the right), where it encounters the 
boundary S. For the £ — mode the potential at the 
boundary is discontinuous and negative, and the field 
now changes sign. This sign change is the large dip in 
Fig. ^ (note that this is a log scale, so that as the field 
passes through ip = the logarithm goes to negative in- 
finity ). It is when this information reaches the event 
horizon that the evolution of the field deviates from its 
evolution in pure Schwarzschild spacetime. The field con- 
tinues to decay with a power-law falloff, but the falloff is 
much slower than in the pure Schwarzschild case. 

The discussion for the £ = 1 case is similar, until the 
reflected wave reaches the boundary E. This is because 
the potential at the boundary is discontinuous and pos- 
itive for ^ > (see equation 10b). Therefore the field 
will be partially transmitted through the barrier at the 
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FIG. 5. Absolute value of the field on the event horizon 
as a function of v, for 1 = 0. The solid line is the evolution 
of the the scalar wave in pure Schwarzschild spacetime. The 
dashed line is the evolution in the Schwarzschild-Einstein-de 
Sitter spacetime. The sharp dip near v = 1000 is where the 
field learns about the boundary E and changes sign. Before 
this point, the field decays as ip ^ . After this point, the 
field decays as ip . 
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FIG. 6. Absolute value of the field on the event horizon, as 
a function of v for £ — 1. The solid line is the evolution in pure 
Schwarzschild spacetime. The dashed line is the evolution 
in the Schwarzschild-Einstein-de Sitter spacetime. The inset 
is a close-up of the region near v — 1000 where the field 
learns about the boundary E. Echoes of the quasi-normal 
oscillations in the field can be seen. 
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FIG. 7. Evolution of the scalar field. Here the boundary 
between the two spacetimes is the line E, with the black hole 
region to the left and the cosmological region to the right. 
The line marked V represents the maximum of the poten- 
tial. The dashed lines are the reflection and transmission 
of the wave pulse from the potential and the boundary. At 
the event horizon the field initially behaves as if it were in 
pure Schwarzschild spacetime (Schild), then evolves differ- 
ently (SFRW). 



boundary and partially reflected off. The transmitted 
wave will make its way off to future null infinity. As for 
the part of the wave packet that has now been reflected 
twice, it will fall back towards the black hole where it 
once again encounters the localized potential. The part 
of the wave that manages to make it through the poten- 
tial on its second encounter heads back towards the event 
horizon, carrying information about the existence of the 
boundary. This second encounter with the localized po- 
tential has the same effect on the packet as it did the 
first time — namely, the field again exhibits quasi-normal 
ringing (see inset of Fig. This "echoing" phenomenon 
occurs only for the ^ > modes of the field. 
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